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ABSTRACT 

Selection of galaxy clusters by mass is now possible due to weak gravitational lensing effects. It is an important 
question then whether this type of selection reduces the projection effects prevalent in optically selected cluster 
samples. We address this question using simulated data, from which we construct synthetic cluster catalogues both 
with Abell's criterion and an aperture-mass estimator sensitive to gravitational tidal effects. The signal-to-noise 
ratio of the latter allows to some degree to control the properties of the cluster sample. For the first time, we apply 
the cluster-detection algorithm proposed by Schneider to large-scale structure simulations. 

We find that selection of clusters through weak gravitational lensing is more reliable in terms of completeness 
and spurious detections. Choosing the signal-to-noise threshold appropriately, the completeness can be increased 
up to 100%, and the fraction of spurious detections can significantly be reduced compared to Abell-selected cluster 
samples. 

We also investigate the accuracy of mass estimates in cluster samples selected by both luminosity and weak- 
lensing effects. We find that mass estimates from gravitational lensing, for which we employ the ^-statistics by 
Kaiser et al., are significantly more accurate than those obtained from galaxy kinematics via the virial theorem. 



1. Introduction 
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^^ Observed samples of galaxy clusters are notoriously contami- 

0^ nated by projection effects (Frenk et al. 1990) even when se- 

^^ lected with respect to their X-ray emission (Bartelmann & Stein- 

Pj metz 1996; van Haarlem, Frenk & White 1997; Cen 1997). This 

O is particularly annoying because clusters constitute an impor- 

I tant class of cosmological objects, and cleanly defined cluster 

O samples could provide a wealth of cosmological information. 

-4— > Not only the definition of clusters within samples, but also ob- 

^ servable cluster properties like their richness. X-ray luminosity, 

.. gas fraction, degree of substructure, velocity dispersion and the 

J> mass estimates derived therefrom can substantially be biased by 

k> projection effects. 

'V^ While not solving the principal problem of observing aspheri- 
d cal three-dimensional objects in projection, selection of clusters 
■ ■ ■ by mass rather than luminosity might prove a substantial step to- 
wards more reliably defined cluster samples. With the advent of 
weak-lensing techniques and wide-field surveys, mass-selected 
cluster samples are beginning to come within reach. The tidal 
field of the total gravitating cluster mass determines the amount 
of distortion observed in images of background galaxies (Kaiser 
& Squires 1993). A variety of elaborate tools has been devel- 
oped in the past few years to detect clusters and derive their 
masses utilising distortions of faint-galaxy samples (Fahlman et 
al. 1994; Kaiser 1995; Kaiser, Squires & Broadhurst 1995; Seitz 
& Schneider 1995; Bartelmann et al. 1996; Seitz & Schneider 
1996; Squires & Kaiser 1996 to name just a few). In particular, 
weighted second-order measures of the distortion field in aper- 
tures like the aperture mass statistics (Map-statistics) provide 
means to derive the signal-to-noise ratio 5 of a measurement 
from the data themselves (Schneider 1996). A possible strategy 
to use aperture measures like Map for cluster detection is to move 
the aperture across a wide data field and monitor Map and S as 
a function of aperture position. Points or regions of high S can 
then be identified as potential clusters, leading to the definition 
of mass-selected cluster samples. This technique led to a cluster 
detection in the ESO Imaging Survey (T. Erben, private com- 
munication). A further spectacular example of galaxy clusters 
detected through weak lensing is given by Kaiser et al. (1998). 



Strictly speaking, distortions caused by gravitational lensing 
trace the gravitational tidal field rather than the mass itself. This 
is due to the fact that a sheet of constant surface mass density 
added to any given lens does not change the shapes of distorted 
images, but only their sizes (the so-called mass-sheet degener- 
acy). Therefore, cluster detection techniques based on distortion 
measurements alone select clusters not by mass, but by the am- 
plitude of their tidal field. Keeping this in mind, we will here- 
after speak of mass-selected clusters in that sense. 

The question then arises, are the cluster samples so obtained 
any more reliable than, e.g., cluster samples selected via Abell's 
criterion or X-ray luminosity? More precisely, what fraction of 
true three-dimensional clusters is detected that way, and what 
fraction of the mass-selected samples are spurious detections, 
either not corresponding to real clusters at all or to clusters out- 
side the desired mass range? This is the question addressed in 
this paper using simulated A^-body data, to which we apply for 
the first time the Map statistics to construct synthetic cluster sam- 
ples. 

We use high-resolution large-scale structure simulations and 
first identify three-dimensional clusters in real space. We then 
populate the simulation volume with galaxies, fixing the aver- 
age mass-to-light ratio. In projection, we apply Abell's criterion 
to construct optically-selected cluster samples for comparison, 
and the S'-statistics to define mass-selected cluster samples via 
the gravitational lensing effects of the simulated matter distribu- 
tion. The quality of the mass-selected samples is then assessed 
in comparison to that of the optically-selected samples. More- 
over, we derive masses both from galaxy kinematics and weak 
lensing to compare the reliability of the different mass estimates. 

Section 2 details the methods used. Section 3 discusses the 
quality of the cluster samples in terms of completeness and the 
fraction of spurious detections. In Sect. 4 and in the Appendix, 
we give examples for the line-of-sight structure of some repre- 
sentative mass-selected clusters. Mass estimates are presented 
and discussed in Sect. 5, and Sect. 6 summarises our conclu- 



Methods 



2.1. N -body simulation 

In order to study the influence of projection effects on cluster 
surveys selected by optical and gravitational lensing informa- 
tion, we need simulated data allowing to mimic as accurately as 
possible the selection of clusters and the determination of their 
properties, for instance their masses. At the same time, the full 
phase-space information is required to assess the amount of con- 
tamination of selected clusters by intervening matter along the 
line-of-sight (hereafter los). For this purpose, we use a large, 
high-resolution A^-body simulation of a standard Cold Dark Mat- 
ter (SCDM) universe. The simulation was carried out within 
the GIF project (Kauffmann et al. 1998) using a parallelised 
version of the HYDRA code. HYDRA (Couchman, Pearce & 
Thomas 1995; Couchman, Thomas & Pearce 1996) is an adap- 
tive particle-particle particle-mesh (AP^M) code. It uses direct 
force summation in clustered regions, whereas long-range forces 
and forces within weakly clustered regions are computed on a 
mesh. The simulation used here transports 256^ particles in a 
periodic comoving cubic volume of (85/z^^Mpc)^, where h is 
the Hubble constant in units of lOOkms^' Mpc^^ 

The simulation adopts the approximation to the linear CDM 
power spectrum (Bond & Efstathiou 1984) given by 



P{k) = 



Ak 



[l + N+(/7<7)3/2+(c^)2]vp/v 



(1) 



where q — T ^ k with the shape parameter F, a = 6.4 h ' Mpc, 
b — 3/z^^Mpc, c = 1.7/z^^Mpc, and v = 1.13. The normalisa- 
tion constant. A, is chosen by fixing Og, the rms density contrast 
in spheres of 8/z^'Mpc radius. It is determined following the 
procedure outlined by White, Efstathiou & Frenk (1993) to meet 
the present-day local cluster abundance of « 8 x 10^^/z^ Mpc^^ 
for rich galaxy clusters. The cosmological parameters are Hq = 
50kms-' Mpc-i (i.e. h = 0.5), i2o = 1, i^A = 0, F = 0.5, and 
08 = 0.6. 

We select a simulation box located at a redshift of z = 0.431 
to achieve high lensing efficiency on sources at redshifts around 
unity, where we assume sources to be throughout this paper. 

For the analysis of gravitational lensing effects of the sim- 
ulated matter distribution, namely the S- and ^-statistics to be 
introduced below, the high resolution provided by the GIF sim- 
ulations is essential. Spatial and mass resolution must be dis- 
tinguished. The spatial resolution is determined by the comov- 
ing force softening length, Zjoft = 36/z^'kpc. At one softening 
length from a particle, the softened force is about half its Newto- 
nian value. This limitation is reduced by the high redshift of our 
simulation box, where the force softening translates to a very 
small force softening angle, 850^ w 10". The mass resolution, 
which describes the effect of the finite particle number, is given 
by the particle mass mp = 1.0 x lO^^/z^^M©. The finite mass 
resolution introduces a white noise component into the simula- 
tions. This is not negligible for the SCDM model because there 
a higher proportion of the particles is in voids than for models 
with lower mean density. 

Since we want to evaluate two different methods for detecting 
clusters in projection, we first have to create a sample of true 
3-dimensional (3-D) clusters or groups from the simulation that 
will serve as a reference set. 

We extract clusters and groups from the 3-D dark-matter sim- 
ulation with a friends-of-friends algorithm (also called group- 
finder, cf. Davis et al. 1985). The friends-of-friends algorithm 
is based on a percolation analysis: It identifies groups and clus- 
ters in the simulation box by linking together all particle pairs 
separated by less than a fraction b of the mean particle separa- 
tion. Each distinct subset of connected particles is then taken 



as a group or cluster. We have chosen b — 0.2, but the result of 
group finding does not sensitively depend on the exact choice of 
b. 

Assuming that the clusters and groups found by the group- 
finder are completely virialised, we compute their virial masses 
M200, defined as the mass enclosed by a sphere with a radius 
r2oo which contains a mean overdensity of 200 times the crit- 
ical mass density, p = 200pcrit- For an ilo == 1 universe, this 
radius approximately separates virialised regions from the infall 
regions of the haloes (Cole & Lacey 1996). Several reference 
sets of "true" 3-D clusters are then formed by selecting objects 
with M200 above certain mass thresholds. 

2.2. Construction of mock cluster catalogues by optical 
cluster selection 

Having dark-matter particles only, we need to populate our sim- 
ulation with galaxies for optical cluster selection. We employ 
the following scheme. 

Galaxy luminosities L are drawn from a Schechter function 
(Schechter 1976), 



(j)(L) = A^* (L/Uy°- exp{-L/U) , 



(2) 



with parameters L* =3.77 x 10 Lq and a = 1 taken from the 
CfA redshift survey by Marzke, Huchra & Geller (1994). The 
formal divergence for L — > in the number-density integral of 
the luminosity function is avoided by introducing a lower lumi- 
nosity cut-off Lo = 0. 1 L* . For the normalisation of the luminos- 
ity distribution, we follow the prescription by Schechter (1976). 
We calculate a richness estimate by computing the most proba- 
ble value of the third-brightest absolute magnitude Mj, and then 
integrate the luminosity distribution from M3 to M3 + 2. Frenk 
et al. (1990) showed that this yields the dimension-less normal- 
isation factor N^, = 60.0. The normalisation factor determines 
the amount of luminous galaxies to be introduced into the simu- 
lation. The total mass-to-light ratio of the 3-D clusters turns out 
to be M/L — 3OO/zM©/L0 on average, in qualitative agreement 
with observations. 

Assuming that mass follows light in our model universe, 
galaxies inherit positions and velocities from randomly selected 
dark-matter particles. In this sense, our constructed galaxy sam- 
ple is unbiased both in number density and velocity. 

Transforming luminosities to apparent magnitudes for higher 
redshifts, we account for the k correction. If the spectral energy 
distribution varies with frequency v as a power law with expo- 
nent a', the additive k correction is 



fc = -2.5(l+a')logio(l+z) 



(3) 



We choose a' = — 1 .5 for the spectral index, which sufficiently 
well reflects the spectral properties of ordinary galaxies. 

Volume-limited cluster catalogues are then obtained after pro- 
jecting particle positions onto planes along the three orthogonal 
axes of the simulation box. Groups and clusters in projection are 
identified with a 2-D version of the friends-of-friends algorithm. 

We then apply the optical Abell criterion (Abell 1958) to se- 
lect galaxy clusters. This widely used cluster detection and clas- 
sification scheme does not depend on redshift z- Therefore, it 
can also be employed at the high redshift of the simulation even 
though it has traditionally been used only for fairly shallow clus- 
ter surveys. Briefly, a cluster is classified as an Abell cluster if 
within the Abell radius of ra = 1 .5 h^ ' Mpc from its centre, and 
after subtraction of the mean background, the number count of 
galaxies exceeds a certain value n^. Counting is restricted to the 
apparent magnitude interval {mj,,mj +2), where mj denotes the 
apparent magnitude of the third-brightest cluster galaxy. The 



actual count n^ is used to assign Abell richness classes ^ . For 
2^ — 0, a cluster has to contain at least n^ — 30 galaxies, while 
5^ = 1 and %_ =2 correspond to tia > 50 and «a > 80, respec- 
tively. 

We also straightforwardly apply Abell's criterion to three- 
dimensional clusters in order to assess the influence of projec- 
tion effects on richness-class estimates. 

For the background subtraction, we follow Frenk et al. (1990). 
In order to estimate the background, i.e. the contamination by 
foreground and background galaxies in the simulation box, we 
assume that the number of galaxies contributing to the contam- 
ination is proportional to the volume projected onto the cluster. 
In our case, we expect 8 background galaxies within a cylin- 
der of volume Vp^oj = 7t'"n/box- Therefore, a cluster with richness 
class 5^ = 1 has to encompass at least 58 galaxies in the appro- 
priate magnitude interval; 8 background galaxies in addition to 
the 50 genuine cluster members. 

Since observed column densities towards galaxy clusters will 
be considerably larger than assumed here, and since the condi- 
tions in realistic observations are less controlled than here, pro- 
jection effects could even be larger in reality. 

2.3. Detection of dark-matter concentrations through weak 
gravitational tensing 

2.3.1. Basic relations 

We briefly review in this section relations from gravitational 
lensing theory important later on. For a derivation cf. Schnei- 
der, Ehlers & Falco (1992). Some remarks on the numerical 
calculation of lensing properties will also be made. 

The dimension-less surface mass density K (also called con- 
vergence) is given by 



k(x) 



E(x) 



tri-diagonal system of linear equations, and back-transforms to 
real space. In contrast to other approaches, the approximation 
is made here by discretising the equations, which can then be 
solved exactly by a subsequent discrete FFT. 

Having determined the deflection potential x|/(x), the local 
properties of the lens, such as the surface mass density K and 
the complex shear Y= ji +iY2, can be expressed in terms of 
second derivatives of \\i{x), 



k{x) 



1 



vV(-?) 



Yi(^) = -2 (¥.11 -¥.22), 

Y2(x) = -\|/.12 , 



(8) 

(9) 
(10) 



where indices / following commas denote partial derivatives 
with respect to jc,. [] 

2.3.2. Aperture mass measures 

Galaxy clusters can be selected solely by their mass using the 
method developed by Schneider (1996). For detecting dark- 
matter concentrations through image distortions of faint back- 
ground galaxies, we define an aperture mass measure m{x) as 



tix) = Jd^yK(y)W{\x-y\), 



(11) 



where the integral extends over a circular aperture |x — y] < R, 
and W{y) is a continuous weight function depending on the 
modulus of y only, vanishing for y > R. We assume W (y) to 
be a compensated filter function, i.e. 



(4) 



/ dyyW{y)^0. 
Jo 



(12) 



with the critical surface mass density 



D, 



4kG DiDi, 



(5) 



Ds, Di, and D^^ are the angular-diameter distances from the ob- 
server to the sources, from the observer to the lens, and from the 
lens to the sources, respectively. The surface mass density k{x) 
is related to the effective deflection potential \(/(x) through the 
Poisson equation 

2k(x), (6) 



vV(^) 



which can be solved for \(/, 



\|/(x) 



dVK(.x')ln(|x- 



(7) 



Boundary conditions have to be specified when solving eq. (I 
numerically. Periodic boundaries are adequate because of the 
periodicity of the simulation volume. 

For numerically computing \(/, the projected particle posi- 
tions are interpolated on a grid of 2048^ cells to maintain the 
high resolution of the A^-body simulation. The resulting surface 
mass density is scaled with the critical surface mass density (q|) 
to find the convergence K. For a numerically stable and effi- 
cient method to convert K to \|/, we use a fast Poisson solver 
(Swarztrauber 1984). The efficiency of this method rests on the 
fast Fourier transform (FFT) leading to an asymptotic operation 
count of o{2NlogN). The algorithm approximates the Lapla- 
cian on a grid, transforms to Fourier space, solves the resulting 



For such filter functions, the aperture mass measure m{x) can be 
expressed in terms of the tangential shear component yt inside a 
circle with radius R (Fahlman et al. 1994; Schneider 1996) 



i{x) = ld^yQi\y\M^-J) 



where Q is related to the filter function W by 

g(x) = -r / dx'x'W(x')-W(x) 
x^ Jo 



(13) 



(14) 



In the above equations, the tangential shear Jtix^y) at the posi- 
tion y relative to position x is given by 



Yt(x;y)=9^[Y(3/)e-2i*] 



(15) 



where (|) is the position angle of y — x. Equation (Oh relates 
the spatially filtered aperture mass m{x) to the observable shear 
field. 

2.3.3. Signal-to-noise ratio 

An estimate for the shear field y, and thus for the aperture mass 
m{x) via equation (pj|), is provided by the distortions of images 
of faint backgroundgalaxies. The complex ellipticity of galaxy 
images, %, is commonly defined in terms of second moments 
of the surface-brightness tensor (Schneider & Seitz 1995). For 

'Notice the signs of tlie sliear components: We follow the sign convention of 
Schneider & Seitz (1995). 



sources with elliptical isophotes with axis ratio r < 1, the modu- 
lus of the source ellipticities is given as \x^^'\ = (1 — r^)/(l+r^), 
and the phase of the x''^' is twice the position angle of the ellipse. 
It has been demonstrated (Schramm & Kaiser 1995; Seitz & 
Schneider 1997) that the ellipticity X of ^ galaxy image is an 
unbiased estimate of twice the local shear in the case of weak 
lensing, |y| <C 1 . If one assumes that the intrinsic orientations of 
the sources are random, 



(x'^' 



0. 



(16) 



with the average taken over an ensemble of sources, then all 
average net image ellipticities (x) 7^ reflect the gravitational 
tidal effects of the intervening mass distribution. We draw the 
source ellipticities from a Gaussian probability distribution 



PsilX^'^)- 



1 



■ exp 
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(-s)|2^ 



(17) 



7ta2[l-exp(-a^2^] \ -x / 

where the width of the distribution is chosen as a^ = 0.3. 
We set the number density of the background sources to n = 
35arcmin^ . 

The complex image ellipticity % can then be calculated in 
terms of the source ellipticity x'^' and the reduced shear g = 
Y(1 — k)"^ by the transformation (Schneider & Seitz 1995) 



,W. 



X 



-2g + gY'> 



1 



■23i{gX 



>)*\ 



(18) 



In analogy to the tangential shear component yt occurring in 
( 13 ), a similar quantity for the image ellipticities can be defined. 
Consider a galaxy image / at a position x,- with a complex im- 
age ellipticity Xi- The tangential ellipticity Xti(^) of this galaxy 
relative to the point x is then given by Schneider (1996) 



Xt,(^) = 9^U 



{X,~X) 



(19) 



where X/ = x,i + ix^ and X = xi + 1x2 are complex representa- 
tions of the vectors x,- = (x,i ,xa) and x ~ (xi ,^2). 

We are now able to estimate the weighted integral (|l3]) as a 
discrete sum over galaxy images. 



m{x)K—Y,Xtiix)Q{\xi-x\) 



2n 



(20) 



The dispersion Gm f x ' 
found by squaring (GO 
leads to 

0^(x) 



of the aperture mass estimate m(x) is 
I and taking the expectation value, which 






Xi- 



(21) 



Finally, the signal-to-noise ratio S at position x can be written as 

mix) V2 1iXtiix)Qi\xi~x\) 



Six) 



a mix) 



[L.Q^ 



11/2 



(22) 



2.3.4. The 'C^-statistics 



So far, the formalism for aperture mass measures and their 
signal-to-noise ratios is independent of the choice for the weight 
function W . Specialising W now, we are led to aperture mea- 
sures with different merits. Two principal choices for the fil- 
ter function have been suggested in the literature. One leads to 



the ^-statistics proposed by Kaiser (1995) and first applied by 
Fahlman et al. (1994). It gives a lower bound to the average sur- 
face mass density K within a circle inside an annulus by measur- 
ing the distortions of background galaxy images inside the annu- 
lus. The ^-statistics will be used in Sect, ^for constraining the 
masses of clusters detected through their 5-statistics. The piece- 
wise constant weight function for the ^-statistics reads (Schnei- 
der 1996) 



W^ix) = 



TtXt 



for < X < xi 
for /? < X < 00 



1 



(23) 



Inserting this weight function into eq. ( [l3| ) yields 

1 r? r2 (fr 



'ni^ix) = ^iri,r2) 



n rj-r^ 



-W) 



(24) 



where r is the distance vector between the point under consid- 
eration and X, and ri and r2 > ri are the bounding radii of an 
annulus around x. It can then be shown that ^(ri,r2) is related 
to the mean convergence ic(ri ,r2) in the annulus by 



C(ri,r2)=K(n)-K(n,'-2), 



(25) 



ic(ri ) being the mean convergence in the circle with radius ri 
around x. In other words, the ^-statistics constrains the average 
convergence in a circular aperture through the tangential shear 
in an annulus surrounding the aperture. Since K(ri,r2) > 0, 
^(ri , ^2) provides a lower bound to the mean surface mass den- 
sity enclosed by ri . 

As mentioned before, it is possible to use the image elliptici- 
ties X of the background galaxies as unbiased estimates of twice 
the tangential shear, Xt ~ 2yt. Therefore, the integral in ( p4| ) can 
be approximated as a discrete sum over galaxy images. 



C(n,.2)-^ 
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L 



Xti 



(26) 
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In this study, we want to obtain a lower bound to the total clus- 
ter masses. For a meaningful application of the ^-statistics, it is 
important to include the complete cluster into the measurement. 
This can be achieved following Bartelmann (1995). If we apply 
the ^-statistics to a nested set of annuli with radii r,-, I < i < n, 
then the ^-statistics for an annulus with r,- < rj reads 



Si,/ ^ ^i ^i,j ■ 



(27) 



where ic,- = K(r,) and ic,; = ic(r, , ry). On the other hand, the mass 
Mij in such an annulus is the product of surface mass density 
times the area, 

(28) 



M:. 



'^ij'^'j y 



where the area of the annulus is 



-nirj 



(29) 



The crucial point is now that the mass contained within a circle 
of radius R is always the sum of the masses contained in annuli 
with outer radii rj < R, irrespective of how the area is decom- 
posed into such rings. Keeping this in mind, eqs. ( pTj ) and (E8h 
can be combined into a system of (n — 1 ) linear equations wirnn 
unknowns m^, where the m^ denote the masses in adjacent rings 



The fact that there is one equation less than the number of un- 
knowns reflects the scaHng invariance of the surface mass den- 
sity K. Assuming that the outermost annulus does not contain 
any significant convergence, i.e. M„_i „ = 0, we finafly arrive 
at the following set of equations for the masses Mj enclosed by 



radii ry: 



Mi 



(1 



Ai ■ 



m\ -AijCij 



where m\ is shorthand for 



m\ 



A\,n—Aij^^\ 



(30) 



(31) 



Of course, a lower bound to the total cluster mass could also 
be obtained by placing an annulus around the entire cluster and 
applying the ^-statistics to that annulus only rather than to a 
set of nested annuli. Our approach has two advantages; first, 
it yields a profile of K which allows to assess the location of the 
outer cluster boundary, for which we found that R = 1 .8/1^ 'Mpc 
is an appropriate choice. Second, it uses galaxy ellipticity mea- 
surements in all annuli rather than the outermost only, thus re- 
ducing the noise. However, the errors in the Mj are correlated 
at successive radii, making an immediate interpretation of the 
significance at any given radius less transparent. 

2.3.5. The S-statistics 

Since the ^-statistics is not designed for detecting mass concen- 
trations, its filter function is not optimised for achieving high 
signal-to-noise ratios, leading to high noise levels in a signal-to- 
noise map. Schneider (1996) solved this problem by introducing 
the smooth, continuous weight function 

1 



Ws(x) = < 



for 
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Q<x<ViR 
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(32) 



In the following, the term S'-statistics refers to the signal-to- 
noise ratio obtained from eq. ( 22 ) using the filter function Ws (x), 
which guarantees low noise in tne signal-to-noise ratio map. The 
parameters a, b, and c are determined once Vi and V2 are spec- 
ified; see Schneider (1996). We choose Vi = 0.05 and V2 = 0.8 
in order to achieve high signal-to-noise ratios and evaluate the 
5-statistics for an aperture size of 2 arc minutes. 

3. Completeness of cluster catalogues 

We are now in a position to investigate completeness and homo- 
geneity of cluster catalogues constructed with the S'-statistics as 
opposed to the optical Abell criterion. To this end, we create 
two different samples of 2-D clusters by applying both meth- 
ods to simulated data projected along the x-, y-, and z-axes. We 
then compare these 2-D clusters with our reference set of 3-D 
clusters. 

To assess the quality of the constructed 2-D catalogues, we 
use several reference sets of 3-D clusters with different mass 
ranges. Looking at different mass ranges instead of cluster rich- 
ness estimates is motivated by two reasons. First, the physical 
quantities of primary interest are the masses. Furthermore, this 
kind of comparison is more suitable for the 5-statistics, which 



does not depend on the distribution of luminous galaxies but 
of the dark matter only. This way of addressing projection ef- 
fects in cluster catalogues differs from previous studies (e.g. Cen 
1997; van Haarlem et al. 1997), which focused on the influ- 
ence of projection effects on the richness estimate of Abell cat- 
alogues. 
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Fig. 1. — Performance of Abell's criterion in identifying clusters in dif- 
ferent mass ranges (cf. Tab. hi). The hatched bars illustrate the number 
of 3-D clusters for each mass range, the solid and dotted bars the num- 
ber of Abell clusters with 3?. > and 3^ > 0, respectively. The overlap 
between hatched and empty bars illustrates the fraction of correct detec- 
tions, and the fraction of spurious detections by Abell's criterion corre- 
sponds to the rest of the empty bars. The fraction of correctly identified 
3-D clusters and the fraction of spurious detections both increase for 
more massive cluster samples. 

The results for the Abell-selected cluster catalogues are sum- 
marised in Tab. M and Fig. m The first column of Tab. |l] lists 
the mass range of the 3-D cluster reference set. The next two 
columns show the fraction of 3-D clusters correctly detected by 
Abell's criterion, and the fraction of 2-D objects which do not 
correspond to 3-D clusters within the chosen mass range, respec- 
tively. The fraction of detected clusters is given with respect to 
the 3-D clusters in the given mass range, while the spurious de- 
tections are given relative to the total number of 2-D detections. 
The last two columns show the same information for a larger 
sample of Abell clusters also including clusters of richness class 
!!{_ =0. Figure |l| illustrates the information of Tab. |l] as a his- 
togram. 

The first mass range considered in Tab. |l| (1.32 — 3.5) x 
IO'^/z^'Mq, reflects the masses expected for Abell clusters with 
richness class 5^ > 0. Looking at absolute numbers, we find that 
the total number of 2-D clusters with ?^ > in the Abell cata- 
logue is very similar to the number of 3-D objects in this mass 
range (in Fig. |l|, the corresponding bars are comparably long). 
However, as Tab. n^ shows, only 53% of the 3-D clusters from 
the reference set can be found in the 2-D sample of i^ > clus- 
ters. On the other hand, a high percentage (50%) of the 2-D 
Abell clusters does not correspond to a true 3-D object. This 
means that not only half of the 3-D clusters in this mass range 
are missed by Abell's criterion, but also a lot of spurious 2-D 



Table 1 . — Completeness and homogeneity of optically selected, synthetic cluster catalogues constructed by Abell's criterion. For 
five mass ranges of 3-D clusters, the detection rate (det.) and the spurious detection rate (spur, det.) is shown for cluster richness 
5^ > and 'J{_ >0. The detection rate gives the percentage of 3-D clusters correctly identified by Abell's criterion. The spurious 
detection rate is the fraction of detected clusters which do not correspond to true 3-D clusters within the considered mass range. 
Further explanation is given in the text. 

3-D clusters ^ >0 ^ > 

Mass range in IO^'^/z^'M© det. spur det. det. spur det. 



1.32- 


-3.50 


53% 


50% 


66% 


82% 


1.03- 


-3.50 


27% 


50% 


56% 


69% 


0.82- 


-3.50 


17% 


50% 


44% 


64% 


0.55- 


-3.50 


13% 


25% 


36% 


40% 


0.10- 


-3.50 


13% 


6% 


32% 


29% 



objects are found. This is due to two competing effects occur- 
ring in projection. Intrinsically rich clusters may disappear in 
the background, while the richness class of poor clusters can be 
enhanced by small groups and field galaxies collected along the 
line of sight. In the above case these two effects approximately 
cancel, so that the total numbers are approximately correct. 

These results are consistent with the findings of Frenk et 
al. (1990) for projection effects in CDM-like universes with dif- 
ferent biasing parameters b. Comparing Tab. 2 in their paper 
with our results for Abell ^ — I clusters, we find similar projec- 
tion effects for our model universe and their CDM-like universes 
with biasing parameters between b = 2.5 and b = 1.6. A direct 
comparison is difficult because of the different normalisations of 
their model universe and ours. Furthermore, the redshift depen- 
dence of the biasing parameter b is not known, further compli- 
cating a detailed comparison. 

For 2-D galaxy clusters with richness class ^ =0, the frac- 
tion of detected 3-D clusters is slightly increased from 53% to 
66%. At the same time, the number of spurious 2-D clusters, 
i.e. clusters which cannot be linked to 3-D objects in the refer- 
ence set, is increased by more than 30%. A detailed analysis 
of the line-of-sight structure of these clusters reveals that this 
large number of spurious detections is partly due to the additive 
projection of poorer groups corresponding to lower-mass 3-D 
objects. This increase in the number of both detections and spu- 
rious detections reflects the enhancement or reduction of cluster 
richness classes due to projection. 

Extending the reference set of 3-D clusters to lower masses 
substantiates this assumption. The number of spurious detec- 
tions declines quite steeply from over 80% for the mass range 
of (1.32 - 3.5) X 10^"^ h-^MQ to below 30% for a lower mass 
threshold of 0.1 x IO'^^/i^'M©, which is nearly one order of 
magnitude smaller than the lower mass threshold for ^ >0 clus- 
ters. Therefore, many of the 2-D clusters detected by Abell's 
criterion do indeed correspond to true 3-D objects, but in very 
different mass ranges. This clearly indicates that for our model 
universe a change of the richness estimate due to projection is 
likely. But still the number of truly spurious detections, i.e. de- 
tections of 2-D objects which cannot be connected with any 3-D 
object, remains quite high even in the broadest mass range. 

Turning to the performance of the 5-statistics in constructing a 
complete and homogeneous catalogue. Tab. ^jand Fig. Q display 
the results of the 5-statistics in a manner analogous to Tab. Iljand 
Fig. [l]for Abell's criterion. Again, the first column contains the 
mass range of the investigated 3-D reference set, while the fol- 
lowing columns display the percentage of detected 3-D clusters 
and of spuriously detected 2-D objects above a certain S'-value. 
The analysis is performed for objects detected above different 
signal-to-noise thresholds. 

The first 5-threshold investigated in detail is S > 5. This 



Mass Range 
(lO'^h-'Mg) 



0.10-3.50 



0.55-3.50 



0.82-3.50 



1.03-3.50 



1.32-3.50 




Fig. 2. — Performance of the 5-statistics in identifying clusters in dif- 
ferent mass ranges (of. Tab. Q). See the caption of Fig. [fl for the meaning 
of the overlapping and non-overlapping parts of the bars. The different 
line types are for different 5 thresholds, as indicated in the plot. With 
increasing 5, the fraction of spurious detections and the completeness 
are both reduced. 



value has been advocated in the Uterature, e.g. by Schneider 
& Kneib (1998), as a signal-to-noise ratio promising signifi- 
cant detections. In comparison to the optical Abell criterion, 
the 5-statistics has a similar detection rate for Abell ^ > like 
objects (53% with Abell's criterion compared to 60% with the 
S'-statistics). The number of spuriously detected objects (50%) 
is identical to that for Abell's criterion in the highest mass range. 
The differences between the two methods show up when de- 
tections and spurious detections at mass ranges with a lower 
mass threshold are considered. Looking at the detection rate 
of spurious 2-D objects, we see a much steeper decline as in the 
Abell case. For a lower mass threshold of 1.03 x lO'^/z^^M©, 
only 27% of the 5-detected clusters do not correspond to 3-D 
clusters of the reference set, whereas more than half of the Abell 
clusters in that mass range are spurious detections. For an even 
lower mass threshold of 0.82 x lO'^/i^ ' Mq, the rate of spurious 
detections falls to only 11%, which clearly indicates that a large 
number of suspected spurious detections in reality corresponds 



Table 2. — Completeness and homogeneity of a catalogue constructed with the 5-statistics based on a weak gravitational lensing 
analysis of distorted images of background sources. Thepercentage of detections (det.) and spurious detections (spur det.) is given 
for different 5-values. For further explanation, see Tab. njand the text. 

3-D clusters S>5 S>4.5 S>4 S>3.5 S>3 



Mass range in 10 /j M© det. spur. det. det. spur det. 



det. 



spur det. det. spur det. det. spur. det. 



1.32-3.50 
1.03-3.50 
0.82-3.50 
0.55-3.50 
0.10-3.50 



60% 


50% 


60% 


68% 


100% 


65% 


100% 


80% 






43% 


27% 


60% 


37% 


76% 


46% 


93% 


63% 






33% 


11% 


48% 


24% 


64% 


32% 


91% 


47% 






18% 


6% 


28% 


10% 


37% 


21% 


61% 


28% 


81% 


32% 


14% 


6% 


22% 


7% 


30% 


16% 


46% 


29% 


90% 


22% 



Mass Range 
(10"h-'M^) 



0.10-J.50 



0.55-3.50 



0.82-3.50 



Rao 



S64 



ZM 



more than one peak coming from within a plateau of high S (see 
Fig. ^. It is important to properly categorise these structures as 
belonging to a single 3-D object and to not count them twice. An 
example of such a situation is shown in Fig. Q. The contour plot 
shows a blow-up of the S-map around two peaks which almost 
overlap in the lower-resolution S-map of the whole simulation 
box (see the mark in Fig. 0). These two peaks in the 5-map 
correspond to one of the most massive clusters in the simula- 
tion. This is reflected in the high 5 = 7.6 of the higher peak, 
while the second peak has S = 3.6. On the other hand, the sam- 
ple also contains examples of 3-D clusters projected onto each 
other showing only a single featureless peak. Therefore, we con- 
clude that the morphological information contained in the S-map 
is low. Both cases of S'-signals will be discussed in more detail 
in Sect. |]. 



1.03-3.50 



1.32-3.50 



Fig. 3. — Comparison of the performances of the 5-statistics and 
Abell's criterion in identifying clusters in different mass ranges. See 
the caption of Fig. Ill for the meaning of the overlapping and non- 
overlapping parts of the bars. The solid and dotted bars are for Abell- 
clusters with 3(, > and 5-selected clusters with 5 > 4, respectively. 
Evidently, the 5-statistics performs better in all mass ranges in terms of 
spurious detections, and completeness is generally also larger. 



to 3-dimensional matter concentrations of lower mass. 

If we reduce the S'-threshold for significant detections to 
eg- S > 4.5 or even below, the detection rate of 3-D clusters 
increases, which means that the S'-cluster catalogue becomes 
more complete. However, the trade-off for the completeness 
is a higher number of spurious detections, belonging to lower- 
mass 3-D mass concentrations. For an 5-threshold of 4.0, we 
are able to construct a catalogue containing all massive Abell- 
like clusters at the expense of also detecting many less massive 
3-D objects. Figure |3] compares the performance of Abell's cri- 
terion with 5^ > and the 5-statistics with S >4. It shows that 
5-selected cluster samples contain fewer spurious detections and 
are generally more complete than Abell-selected samples. 

A crucial point in the application of the 5-statistics is the iden- 
tification of peaks in the 5-map. Following Schneider (1996), 
we use a circular aperture for the 5-statistics, which leads to an 
increased sensitivity for round objects. However, some of the 
5-maps for our simulation data show extended, non-circular ar- 
eas with significant 5-signals. Several of these structures contain 




40.4 41.2 42.0 42.8 
X in Mpc/h 



43.6 



Fig. 4. — 5-map for the rectangular section in Fig. pi showing the 
double-peak structure of the map in more detail. Contours are spaced 
by A5 = 1. The structure features two maxima, but corresponds to a 
single 3-D cluster. 

Summarising the quantitative results from both cluster detec- 
tion methods, we can say that the 5-statistics leads to better re- 
sults than Abell's criterion. The catalogues constructed with the 
5-statistics are more complete and suffer less from spurious de- 
tections, at least in the sense that most peaks coiTespond to true 
3-D objects. The 5-statistics evidently produces fewer truly spu- 
rious detections than Abell's criterion. However, we note that it 
is not possible to obtain a complete catalogue by counting only 
peaks with a high signal-to-noise value 5 > 5. There is no strict 



correlation between the height of a signal in the 5-statistics and 
the mass of identified 3-dimensional objects. 



u 



C 
>- 



80 


- 


' 1 


* • 




' ■ ■ fc' _ 


60 


9 


« 


o 

0° o 


ffiO 


.ft ■ 

■ 




40 




♦o 


m 


# 


• 


□ 

o 


o 

0* 




■ ^ 




m 




■ 


20 


- 








•• 








o 












i^ 











1 


1 




1 .n. , 1 



20 40 60 80 

X in Mpc/h 



Fig. 5. — Projection of true 3-D clusters/groups (squares) along the z- 
axis onto the full xy-plane of the simulation box, and comparison with 
Abell clusters identified in projection (2-D-Abell clusters, circles), and 
mass-selected clusters identified with the 5-statistics (diamonds). 3- 
D clusters with masses in the range of (0.10-3.5) x 10^'^ h^^ Mq, 2- 
D clusters with ^ > (Abell), and 5 > 3.0 (Map) are shown. Filled 
symbols refer to clusters with masses > 0.55 X 10 /i^'Mq,^ >Oand 
5 > 3.5. (The offset between close symbols is deliberate to facilitate 
reading; it does not reflect the accuracy of measuring cluster positions.) 



4. Structure of representative clusters 

To achieve a deeper understanding of projection effects, we 
study the structure of archetypical clusters or groups along the 
line-of-sight. The main emphasis in this section will be on clus- 
ters selected by the S'-statistics. A detailed discussion of the 
structure of Abell selected clusters along the line-of-sight for 
both real space and velocity space can be found, e.g., in Cen 
(1997) and van Haarlem et al. (1997). 

Even though this section will concentrate on clusters detected 
by the 5-statistics, in some cases also Abell-selected clusters will 
be discussed if the S'-selected clusters also satisfy Abell's crite- 
rion. As Fig. shows for one of the three projection directions, 
this is the case for a lot of S'-selected objects, i.e. there is consid- 
erable overlap between the two selection methods. Abell's crite- 
rion detects the visible light from galaxies, while the 5-statistics 
is sensitive to the underlying distribution of dark matter, mak- 
ing it possible to construct a "mass-selected" sample of clusters, 
as opposed to "flux-limited" samples which are obtained by ob- 
serving luminous galaxies. Since this study is performed un- 
der the supposition that both selection methods detect the same 
physical objects, we have to assume that luminous galaxies are 
good tracers of the dark matter distribution. This is secured by 
the assumption used to populate the simulation that light follows 
mass. 

For a more detailed analysis, the 5-selected clusters will be 
subdivided into three classes according to the 5-threshold em- 
ployed. The first class considered are clusters detected with 
S > 5, the next class contains clusters with 5 > 5 > 4, and the 



last class clusters or groups with 4 > 5 > 3. The division into 
classes according to signal-to-noise values allows an investiga- 
tion of systematic differences of projection effects in the differ- 
ent classes. 

It is expected from theory that higher-mass clusters lead to 
larger values in the 5-map. Such a trend was found in Sect. b1 
but there is no sharp correlation between the masses of detected 
3-D clusters and the threshold imposed on the 5-statistics. This 
can be explained in terms of the intervening matter along the 
line-of-sight. Lower-mass 3-D objects are more prone to pro- 
jection effects in the sense that the intervening matter makes up 
a more substantial fraction of their mass. Therefore, projection 
effects become more important for clusters with lower S. In the 
case of clusters with intrinsically lower masses, less intervening 
matter is needed to alter the signal in the 5-map. 

4.1. S-statistics: S >5 

As discussed in Sect. ^ there is a good correspondence between 
5 > 5 clusters and massive 3-D-clusters. Investigating the line- 
of-sight structure of these 2-D-clusters, we can generally state 
that nearly all of them show a high, pronounced peak in the po- 
sition histogram at the position of the true 3-D cluster. Even 
though the amount of contamination with intervening matter in 
this group is only moderate, some velocity histograms deviate 
significantly from a Gaussian shape. 

The cluster given in Sect. |] as an example for a cluster ex- 
hibiting a double-peaked 5-map (see Fig. Q) clearly belongs into 
this class, since the main peak has S — 7.3, while the nearby 
second peak has only 5 = 3.6. The structure along the line-of- 
sight in both real space and velocity space is shown in Fig. ^ 
The position histogram is characterised by a dominant peak 
at the position of the corresponding 3-D cluster which has a 
mass of M3_D = 3.4 x IO'^/j^'M©. The amount of dark matter 
along the line-of-sight is moderate with only two small clumps 
40 Z;^ ' Mpc and 65 Z;^ ' Mpc behind the main clump, which both 
have masses smaller than IO'^/z^'Mq. Although the 3-D clus- 
ter has Abell richness class ^ —Q, the projected cluster satis- 
fies the 2-D Abell criterion for 5^ = 1, indicating an inflation 
of richness class. The velocity dispersion 03_d = 775 km s^' of 
the 3-D cluster is reduced to O2-D = 676 km s^' in projection, 
hinting at an asymmetric velocity ellipsoid of the 3-D cluster. 

Although the cluster is only moderately affected by projec- 
tion effects, the los velocity histogram strongly deviates from a 
Gaussian, which is also true for the velocity distribution of the 
3-D-cluster alone. The deviation from Gaussianity, as measured 
by higher order moments of the distribution like the skewness 
S and the curtosis 3C, is 53-d = —0.62 and 3C3-D = 0.47 for 
the 3-D-cluster as opposed to 52-d — 1-57 and 3C2-D = 3.31 in 
projection. The substantial increase of skewness and curtosis in 
projection emphasises the influence of projection effects in ve- 
locity space. Together with the decrease of the velocity disper- 
sion and the increase of richness class in projection, this hints at 
the presence of non-virialised sub-clumps in the vicinity of the 
main cluster. Yet, this detection coiTesponds to a very massive 
3-D cluster. More examples are given in Figs. O and O] in the 
Appendix. 

4.2. S-statistics: 5 > 5 > 4 

The next sample of 5-detected clusters has lower signal-to-noise; 
but, as has been shown in Sect. ^ the mass-selected cluster cat- 
alogue becomes more complete if these clusters are included. 
Most of them correspond to true 3-D clusters. Even though 
these detections are significant and contain some massive clus- 
ters, their amount of contamination in relation to their main 3-D 
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Fig. 6. — Structure of a large 3-D cluster with moderate contamina- 
tion, detected as a double peak witli S = 7.3 and S = 3.6 (cf. Fig. m. 
It is also an Abell it = 1 cluster, both in projection and in 3-D. Tne 
upper left panel displays the histogram of the dark matter distribution 
in real space along the line-of-sight, and the upper right panel shows 
the corresponding velocity histogram. The two lower panels show the 
dark matter distribution along a cylinder of radius 1.5/?^' Mpc in real 
space (along the line-of-sight) and in velocity space. For better display 
in both lower panels, a random fraction of the dark matter particles is 
omitted. 



cluster is generally larger than for clusters detected with larger 
S. One typical example of this class is shown in Fig. 0, two more 
are given in Figs. [14 and Ol in the Appendix. 
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Fig. 7. — Structure of a moderately contaminated 3-D cluster with S — 
4.7. See the caption of Fig. n for a description of the panels. 

The cluster in Fig. R is a high-mass cluster with M3_d — 
1.5 X 10^^ h^^ Mq, a clump nearby, and a second mass clump 
30/!^' Mpc away. It is detected at S = 4.7. The velocity 
dispersion of the projected cluster is broadened from 03-D = 



697 km s^ to02-D = 865kms^ , and the projected velocity dis- 
tribution has a bimodal shape. The higher-order moments of the 
velocity distribution indicate this through skewness and curtosis 
in projection (52-d = 1-85 and ?(2-d — 3.27), compared to 3-D 
(53-D =0.18 and 3C3-D = -0.40). The cluster satisfies Abell's 
criterion with richness class i**;, = 0, both in projection and in 
3-D. 



4.3. S-statistics: 4 > 5 > 3 

The last class considered contains clusters with S between 3 
and 4. This class is mainly discussed for reasons of complete- 
ness. Clusters identified with such signal-to-noise correspond 
to lower-mass objects making projection effects along the line- 
of-sight more important. This can be seen looking at the three 
examples in Fig. Band in Figs. HQ and O. 
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Fig. 8. — Line-of-sight structure of a highly contaminated cluster This 
cluster is detected at 5 = 3.2. See the caption of Fig. n for a description 
of the panels. 

The peak with 5 = 3.2 corresponding to the cluster shown 
in Fig. g| is due to a 3-D cluster with a rather low mass only, 
M3_D = 5.9 X lO'^/i~^M0. Compared to the low 3-D cluster 
mass, there is a high amount of intervening matter with several 
mass clumps along the line-of-sight. Because of this substan- 
tial contamination, the projected velocity dispersion is overesti- 
mated; 02-D = 669 kms^' compared to 03-0 — 364 kms^^ for 
the 3-D cluster The velocity histogram has a second peak at the 
high-velocity tail of the distribution. The higher-order moments 
of the velocity distribution of the 3-D cluster (53_d = —0.05 and 
aC3-D = —0.07) change by a large amount when looking at the 
projected velocity distribution (52-D = 1-12 and 3C2-D — 2.43). 
This reflects the large influence the intervening matter exerts on 
observation. 

5. Mass ESTIMATES 

The previous two sections put emphasis on the completeness 
of catalogues constructed with two different selection criteria, 
Abell's criterion and the 5-statistics based on weak gravitational 
lensing. Furthermore, we investigated the structure of some de- 
tected clusters along the line-of-sight. Of course, both kinds of 
information are important when deriving statistical information 
from such catalogues. A third very important test for cosmolog- 
ical theories are the different mass estimates and their relation- 
ship with each other A mass estimate closely related to the opti- 
cal selection derives from the virial theorem. As a gravitational- 
lensing based mass estimate, we choose the ^-statistics, w hich is 



closely related to the 5-statistics as demonstrated in Sect. 2.3.4, 
Under the assumption that clusters of galaxies are bound self- 
gravitating systems in dynamical equilibrium, the total clus- 
ter mass can be estimated via the virial theorem (Binney & 
Tremaine 1987; Sarazin 1986), 



Mu 



Rg{v^ 



(33) 



where Rq is the gravitational radius of the cluster relating the 
system's mass to its potential energy. Gunn & Gott (1972) 
showed that this radius is approximately given by r2oo, the ra- 
dius of a sphere containing an overdensity of 200pcrit. For the 
clusters of our study, a radius of Rg — 0.75/i^'Mpc is a good 
approximation. Observationally, only the los velocity dispersion 
Oil can be measured. Assuming isotropic orbits, the two quanti- 
ties are related by (v^) — 30n . 

When calculating the radial velocities from simulated data, 
the Hubble flow has to be added to the peculiar velocities of 
the simulation. Since the simulation data are at high redshift, 
the dependence on redshift of the cosmological constants also 
has to be taken into account. Therefore, the radial velocities are 
given by 



V|| — fl(z)x|| +fl(z)i|| = a{t) [i/(z)x|| +i||] 
with the expansion factor 



and the Hubble parameter 



H{z)=Ho{l+z) 



3/2 



(34) 



(35) 



(36) 



Apart from the validity of the virial theorem, the virial mass 
estimate depends solely on a correct estimate of the velocity dis- 
persion 0|| . Since the velocity dispersion is very sensitive to field 
galaxies and small sub-clumps projected onto the main clusters, 
it is important to remove these from the sample. We convolve the 
velocity histogram with a 4000km s^' wide top-hat filter to re- 
ject interlopers, i.e. galaxies with relative velocities greater than 
4000km s^' from the peak of the convolved histogram are re- 
moved. We further employ the so-called 3o-clipping procedure 
proposed by Yahil & Vidal (1977) which has widely been ap- 
plied to observational samples. It can be summarised as follows: 

1. compute the mean radial velocity; 

2. remove the galaxy which deviates most from the mean of 
the sample and re-determine the mean without this galaxy; 

3. if the removed galaxy deviates from the new mean by more 
than 3o, it is removed from the sample; 

4. repeat the procedure until the last tested galaxy remains in 
the sample. 

Figure ^ displays the correlation of the virial mass estimate 
with the true mass of the corresponding 3-D cluster. The left 
and right panels show the ratio of the virial mass with the 3-D 
cluster mass as a function of the 3-D mass before and after 3a 
clipping, respectively. 

The first thing to notice is that the masses of clusters with rich- 
ness class ?i^ > are less severely overestimated than masses of 
clusters with lower richness. This holds true for the mass es- 
timates before and after 3o-clipping for both the mean and the 
medianjas can be seen in Tab. pt The second thing readily seen 
in Fig. ^ and Tab. g| is the large dispersion of the underlying 
distribution. This dispersion is smaller for clusters with higher 
richness class than for clusters with the lowest richness class 
5^ = considered. We also note that this dispersion is hardly af- 
fected by the 3o-clipping procedure. The only effect of the clip- 
ping procedure is to reduce the average of the estimated cluster 
masses irrespective of the richness class. A third trend to be seen 
in Fig. ^ is that the overestimation of the mass is generally more 
severe for 3-D clusters with lower mass. For the most massive 
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Fig. 9. — Comparison of the virial mass estimate Myj to the true clus- 
ter mass M^-D- Both panels of this figure display the ratio Mvt/^3-d 
as a function of the true cluster mass M^-u- In the left panel, the ve- 
locity dispersions before 3CT clipping are used to compute the virial- 
theorem based mass estimate, while in the right panel, the velocity dis- 
persions after 3o clipping are used. In both panels, the open circles 
refer to 3?_ = clusters and the filled circles to 3J, = 1 clusters. 



clusters in the sample (M3-D < 2.0 x lO'^/i^^M©), the contin- 
uation of this trend in some cases leads to an underestimation of 
the masses, as can be seen in the right-hand side of each panel 
in Fig. ^. The 3o clipping procedure fails to correct for the mass 
overestimates. When comparing Fig. |and Tab. H to Fig. 15 of 
Cen (1997), one has to keep in mind the different selection pro- 
cedure for clusters or groups of galaxies in both studies, but on 
the whole the results are consistent. 

The behaviour described above can largely be attributed to the 
influence of projection effects on the velocity dispersion. Gener- 
ally, the inclusion of field galaxies and unvirialised sub-clumps 
broadens the distribution and leads to distributions which devi- 
ate significantly from Gaussian shape, as illustrated by the ex- 
amples in Sect, ni The clipping procedure is successful when 
the amount of contamination along the line-of-sight is low or 
moderate, but the algorithm fails to remove larger sub-clumps 
projected onto the main cluster which can significantly broaden 
the distribution, sometimes even making it bimodal. In some 
cases it is possible that the clipping procedure removes galaxies 
belonging to the 3-D cluster, thus contributing to an underesti- 
mation of the mass. 

Even though we expect from the studies of Frenk et al. (1990) 
and van Haarlem et al. (1997) that the high-velocity tail of the 
velocity distribution is severely overestimated, the effect on the 
mass estimate is most pronounced for galaxy clusters with lower 
mass. This is due to the fact that they are more easily overesti- 
mated with respect to their true dispersion. 

The ^-statistics as compared to Abell's criterion leads to 
smaller overestimates of the 3-D cluster masses as shown in 
Interpreting the quantitative results of the ^-statistics 



Fig. 10 



mass estimate, one has to keep in mind two competing effects: 
On the one hand, the ^-statistics, like every gravitational-lensing 
based method, measures all the mass along the line-of-sight to 
the cluster; on the other hand, it gives a lower bound to the clus- 
ter mass. In combination, these two competing effects lead to 
fairly moderate mass overestimates, as can be seen in Tab. |} 
This also explains the difference to the lensing mass estimates 
given in the paper by Cen (1997). There, all masses along the 
line of sight are added up under the assumption of a perfect lens- 
ing reconstruction method with an otherwise calibrated mass- 



mean median standard deviation 



Table 3. — Comparison of statistical parameters for the ratio Mvt/A^3-d and Mr/M^^o for different sub-samples of masses 
estimated via the velocity dispersion and the ^-statistics. 
sample 

complete optical sample before 3a clipping 
Abell cluster 'J{_ = 1 before 3o clipping 
Abell cluster ?^ = before 3o clipping 

complete optical sample after 3o clipping 
Abell cluster 1{_ = 1 after 3o clipping 
Abell cluster 1{_ —Q after 3o clipping 

complete lensing sample 
S>5 

5>5>4 
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to the large underlying dispersion in the virial mass estimate. 



3.0 


1 1 1 1 1 1 1 1 


1 


1 1 . 1 . 




''''''''''''''' 


2.5 


<> 




o 




- 


©2.0 


o 








^ 


S 1.5 
9, 

> 




o 
o 




?.* 


o 


a 


> 1.0 


- 


o 


o o 




B 


0.5 


- 




# 




o 

O ( 


0.0 


. . . . 1 . . . 


1 


, , , 1 , 




1 . . . . 1 < < < . 1 . . . . 



0.0 0.5 1.0 1.5 2.0 

M3„ (lO^^h"^ Mq) 



2.5 



3.0 



Fig. 10. — Comparison of the ^-statistics mass estimate to the true 
cluster mass M3_d- Diamonds and circles refer to clusters identified 
with a 5 > 5 and 5 > 4, respectively. Clusters detected at lower S are 
excluded because of their high contamination. 



sheet degeneracy. 

The other interesting feature in Fig. 11^ and Tab. is the low 
dispersion of the underlying distribution. This dispersion does 
not depend sensitively on the S'-value at which the clusters are 
detected. (Clusters detected with 5 < 4 where excluded here be- 
cause of their large contamination.) The dispersion is typically 
less than a third of the dispersion in the Abell samples. 

As for clusters detected with Abell's criterion, masses of 
small 3-D clusters are more strongly overestimated than for 
more massive 3-D clusters. This is due to the fact that the 
proportion of contaminating matter to the 3-D cluster mass is 
higher for less massive 3-D objects than for the extremely mas- 
sive objects. For the intermediate-mass objects, the fact that the 
^-statistics only gives lower bounds to the masses partially out- 
weighs this effect. There, the masses are even slightly underes- 
timated. 

Investigating the relationship between velocity-based mass 
estimates and the gravitational lensing based ^-statistics in 
Fig. O, we see that the ^-statistics gives on average smaller es- 
timates of the 3-D cluster masses than the virial theorem. Again 
we stress that this is due to the fact that the ^-statistics is derived 
under the assumption of an empty outer annulus, restricting it to 
estimate lower bounds to the masses. The dispersion between 
the ratio of ^-statistics mass to virial mass is large, which is due 
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Fig. 11. — Comparison of the ^-statistics mass estimate to the virial 
mass estimate A/yx- Only such clusters are included which are detected 
with both methods, Abell's criterion and 5-statistics. Only clusters with 
5 > 4 are considered. 



6. Conclusions 

We investigated for the first time with simulated data whether 
mass-selected galaxy cluster samples are more reliable than 
samples constructed via Abell's criterion. Selection of clusters 
by mass is possible through their gravitational lensing effects, 
in particular through the coherent image distortion patterns im- 
posed on faint galaxies in their background. As mentioned in 
the introduction, image distortions trace the gravitational tidal 
field of a lens rather than its mass, and it is in that sense that 
we speak of "mass-selected" cluster samples. The second-order 
aperture-mass statistics Map was used which is particularly suit- 
able for detecting dark-matter haloes. The signal-to-noise ratio, 
S, of Map provides a straightforward detection criterion. We also 
compared the performance of cluster mass estimators based on 
cluster-galaxy kinematics and gravitational lensing. Our results 
can be summarised as follows. 

1. As already found in previous studies, Abell clusters are 
severely affected by projection effects. This not only con- 
cerns the selection of Abell clusters, but also mass esti- 
mates based on galaxy kinematics and the virial theorem. 



indicating that the velocity dispersion is also hampered by 
projection effects. A second reason for the failure is the fact 
that the assumption of dynamical equilibrium is not justi- 
fied in at least some of the clusters. The projection effects 
are worse for clusters and groups of lower richness class. 

2. Clusters detected with a high significance S of Map are less 
affected by projection effects than typical Abell-selected 
clusters. Like Abell cluster samples, the mass-selected 
cluster samples are generally incomplete: Samples of clus- 
ters detected above a certain S threshold typically do not en- 
compass all three-dimensional clusters present in the sim- 
ulation; some clusters have lower S. However, the com- 
pleteness of the samples can be increased by lowering the 
S threshold. We therefore investigated the effect of vary- 
ing the S threshold on the samples. Completeness of w 
100% can be achieved for massive three-dimensional clus- 
ter samples (M > IO'^^/i^'Mq) by varying S>4. Then, 
the samples also contain a substantial fraction of spuri- 
ous detections, most of which correspond to real clusters 
with smaller masses. Generally, there is a trade-off be- 
tween completeness and the contamination by spurious de- 
tections. More complete cluster samples are more heavily 
contaminated by spurious clusters, and the balance can be 
adapted choosing the S threshold. It should be noted that 
the exact thresholds on S depend somewhat on the choice 
of the weight function entering the definition of S (cf. the 
discussion in Sect. 2). 

3. While qualitatively the same trend is also observed in 
Abell-selected cluster samples, the S'-statistics generally 
performs significantly better than Abell's criterion: Higher 
completeness can typically be achieved with a lower frac- 
tion of spurious detections. For instance, cluster samples 
detected at 5 > 4 contain all of the most massive clusters 
in the simulation and 65% spurious detections, while Abell 
samples with richness ^ >0 encompass only about two- 
thirds of the most massive clusters and 82% spurious de- 
tections. 

4. Lensing-based mass estimates are significantly more accu- 
rate than mass estimates based on cluster-galaxy kinemat- 
ics and the virial theorem. Virial masses are typically bi- 
ased high because line-of-sight velocity distributions are 
broadened by projection effects. Lensing also adds up 
mass in front of and behind the clusters, but the bias is 
less severe. The standard deviation from the true (three- 
dimensional) mass of the lensing mass estimate is smaller 
by a factor of three or more than that of the virial mass es- 
timates. It should, however, be noticed that the accuracy of 
lensing-based mass estimates depends on the depth of the 
background-galaxy sample and other observational effects. 
While the mass estimates based on the ^ statistics are accu- 
rate to within « 30% in our simulations, they may well be 
less accurate under realistic observational conditions. 

Our study underestimates projection effects because of the 
limited size of the simulation volume. This affects both the op- 
tical and the lensing-based cluster selection. Yet it appears that 
selection of clusters by mass yields more reliable cluster sam- 
ples than optical cluster selection, and, more importantly, the 
quality of the samples can be controlled by an objective crite- 
rion, namely the signal-to-noise threshold imposed. 

Selection of clusters by their gravitational-lensing effects is 
comparable to selection by their X-ray luminosity. In essence, 
both methods measure weighted projections of the Newtonian 
cluster potential along the line-of-sight. However, lensing-based 
cluster detections only require sufficiently deep imaging of wide 



fields in optical or near-infrared wave bands, and detection al- 
gorithms can then be applied in a straightforward manner. In 
particular, lensing can detect clusters out to significantly higher 
redshifts than X-ray surveys. What is more, lensing-based mass 
estimates do not rely on any assumptions on the composition and 
physical state of the cluster matter, in contrast to X-ray mass es- 
timates. It can therefore be expected that reliable, mass-selected 
cluster samples at moderate to high redshifts can be constructed 
in the near future from upcoming deep, wide-field surveys with 
a straightforward, well-controlled algorithm, and that the accu- 
racy of cluster mass estimates will generally be substantially im- 
proved. 

Lensing-based cluster detection algorithms like that based on 
the S'-statistics can be utilised for cosmology without reference 
to actual 3-D clusters. Instead, one could define a "cluster" op- 
erationally as something visible as a sufficiently high peak in an 
S map, and then compare model predictions with observations 
(cf. Kruse & Schneider 1999). If, however, emphasis is laid on 
clusters as individual entities, it needs to be clarified how well 
different selection criteria fare in detecting true clusters. Our 
study has shown that selection of clusters by means of gravita- 
tional lensing techniques can be adapted such that the resulting 
samples are superior to Abell-selected samples in terms of com- 
pleteness, spurious detections, and the quality of mass estimates. 
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A. Structure of further representative clusters 

We give a few more examples of line-of-sight structures of S- 
selected galaxy clusters here. 



A.l. S-statistics: S>5 

A second example for a cluster with high S is given in Fig. |l2[ 
The cluster is detected at 5 = 9.6. The particle distribution in 
real space is broad and dominated by a massive 3-D cluster with 
amass of2.1 x lO^^/j^'Mp^. This cluster is detected as an Abell 
cluster in projection, but the main 3-D cluster by itself already 
passes the luminosity threshold of a 3-D Abell cluster. In con- 
trast to the first example, the velocity dispersion is hardly af- 
fected by projection. The 3-D cluster has a velocity dispersion of 
03_D = 884 kms^', while the dispersion of the projected clus- 
ter is 02-D = 779kms^'. The higher-order moments indicate 
a velocity distribution close to Gaussian shape for both the 3-D 
cluster (i3_D — 0.01, 3C3-D = —0.04) and the projected cluster 
{Si-D = 0.09, %2-D = -0.37). All this reveals a fairly relaxed 
cluster with low contamination. 

Almost all other clusters in this class show similar position 
and velocity histograms. The only exceptions are the 2-D clus- 
ters corresponding to less massive 3-D clusters. For one of these 
clusters with relatively high S = 7.9, the structure is given in 
Fig. O. Even though the position histogram is dominated by a 
3-D cluster, the distribution for this cluster is broad, and there is 
a large amount of intervening matter with at least four smaller 
clumps with masses of order M ~ IO'^/i^'Mq. Qualitatively, 
the los velocity histogram looks artificially broadened by these 
clumps, and in fact the velocity dispersion (as-o = 651 km s^') 
is significantly increased in projection (02-D = 860kms^'). 
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Fig. 12. — Structure of a massive 3-D cluster detected with S - 
See the caption of Fig. a for a description of the panels. 



9.6. Fig. 14. — Structure of a moderately contaminated 3-D cluster with a 
S = 4.7. See the caption of Fig. fe for a description of the panels. 



The higher-order moments are also strongly affected by this in- 
tervening matter (53-d — —0.05 and J<3-d — —0.64 compared 
to 52-D = 119 and ?C2-d = 0.87). This cluster is detected as 
^ = I Abell cluster although it corresponds only to a 3(, = 
cluster in 3-D. 
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Fig. 13. — Structure of a less massive 3-D cluster whose size is in- 
creased in projection due to matter concentrations along the line-of- 
sight. See the caption of Fig. g for a description of the panels. 



A.2. S-statistics: 5 > 5 > 4 



Figure 14 shows a cluster with S = 4.7. It is apparently 
only mildly contaminated by a clump 30/;^'Mpc from the 
main clump, which is a high-mass object with M3_d = 2.3 x 
lO'^/i^^Mpj. The projected velocity dispersion is almost un- 
affected (02-D = 650kms^' compared to 03_d ~ 635kms^'), 
and shows a bimodal feature, which is also reflected by the cur- 
tosis of the projected cluster, 3C2_d ~ 0.51, while the velocity 
distribution of the 3-D cluster has a negative curtosis of ?(3_d — 
—0.22. Similarly, the skewness changes from 52-d = —0.44 to 
53-D = —0.03. The cluster is detected as a 2-D Abell cluster 
with richness class 5(, = 1, while the richness class of the 3-D 
cluster is 5(, = 0. Therefore, the richness class is inflated due 
to projection. Even though this cluster shows some projection 



effects, the corresponding 3-D cluster is massive and therefore 
clusters like that should be included in a mass-limited sample. 
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Fig. 15. — Structure of a moderately contaminated 3-D cluster with 
S = 4.4. See the caption of Fig. W for a description of the panels. 

The last example for this class is shown in Fig. |l5[ Here, 
the 5-map has a peak with S = 4.4. The position histogram 
shows a very broad peak with a secondary maximum on top of 
the main peak. The corresponding 3-D cluster has a high mass, 
M3-D = 3.4 X IO'^^/i^'M©. The projected velocity distribution 
is only moderately skewed with 52-d = 0.22 compared to the 
skewness of the main cluster alone, 53-d = 0.17. However, the 
curtosis of the projected peak, JQ-d = 0.38, even changes sign 
when compared to the 3-D cluster, 9(j-d = —0.27. This cluster 
satisfies Abell's criterion in projection, but the main peak has a 
lower richness class, 5^ = 0. 



A.3. S-statistics: 4 > 5 > 3 

Another example for a low-S' cluster detected at 5 = 3.4 is dis- 
played in Fig. O. This detection also corresponds to a 3-D clus- 
ter with M3-D =6.6 x 1O'^/!"'M0. Again, the velocity distribu- 
tion of this cluster is largely altered by the considerable amount 
of intervening matter. The velocity dispersion itself is inflated 
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Fig. 16. — Structure of moderately large 3-D group with 5 between 3 
and 4. The 3-D object is contaminated by projection along the line-of- 
sight, leading to an increased S of 3.4. See the caption of Fig. |g for a 
description of the panels. 



from a3_D — 504kms^' to O2-D = 1134kms^ . This is re- 
flected by the curtosis, which changes from aC3_D = —0.02 to 
%2-D = —0.5, while the skewness changes from 53-D = 0.18 
to 52-D = —0.94. Both low-S' examples are neither 2-D Abell 
clusters nor do they pass the selection criteria for Abell clusters 
in 3-D. 
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Fig. 17. — Structure of a spuriously detected object which does not 
correspond to a 3-D cluster. The high contamination along the line-of- 
sight leads to 5 = 3.8. See the caption of Fig. o for a description of the 
panels. 



The last example in Fig. |lj with 5 = 3.75 does not corre- 
spond to a 3-D cluster or group with mass exceeding M3-D = 
1O'^/i^'M0. Instead, one sees a large amount of contaminat- 
ing matter and smaller sub-clumps. This material is responsible 
for the signal in the S map. The velocity distribution is char- 
acterised by three peaks with dispersion O2-D = 1235kms^\ 
skewness 52-d = —0.30, and curtosis 5Q-d = —1.39. Obvi- 
ously, the contamination along the line-of-sight is large enough 
to lead to the detection of an Abell cluster with richness class 
^=0. 



